Asymmetric noise probed with a Josephson junction 



o 
o 

(N 



CD 



o 



ON 

in 

o 

o\ 
o 



X 



Q. Le Masne,^ H. Pothier/'B Norman O. Birge,^ C. Urbina,^ and D. Esteve^ 

^ Quantronics group, Service de Physique de I'Etat Condense (CNRS URA 2464), CEA-Saclay, 91191 Gif-sur-Yvette, France 
"^Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 

(Dated: January 16, 2009) 

Fluctuations of the current through a tunnel junction are measured using a Josephson junction. 
The current noise adds to the bias current of the Josephson junction and affects its switching out of 
the supercurrent branch. The experiment is carried out in a regime where switching is determined by 
thermal activation. The variance of the noise results in an elevated effective temperature, whereas 
the third cumulant, related to its asymmetric character, leads to a difference in the switching rates 
observed for opposite signs of the current through the tunnel junction. Measurements are compared 
quantitatively with recent theoretical predictions. 
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The current through voltage-biased electrical conduc- 
tors exhibits fluctuations, which, in contrast to equilib- 
rium Johnson-Nyquist noise, are not symmetric with re- 
spect to the average current. This translates into finite 
odd cumulants in the distribution of the number of elec- 
trons transfered through the conductor in a given time. 
Whereas the full counting statistics of this number can 
be calculated for arbitrary conductors P, [2| , up to now 
high order cumulants have been measured in very few 
experiments. The third cumulant has been successfully 
accessed by signal processing the time-dependent current 
d, d, Q , but with set-ups that are restricted either to low 
impedance samples, which leads to large environmental 
effects or to low frequencies d, @. Another strat- 
egy, put forward by Tobiska and Nazarov consists 
in using a Josephson junctionfJJ) as a large bandwidth 
on-chip noise detector 0, H, |^. It has a high intrinsic 
sensitivity, and can be coupled to noise sources over a 
large range of impedances. The detection principle relies 
on the exponential current sensitivity of the switching of 
a JJ from a metastable zero- voltage branch to a dissipa- 
tive one. When biased at a current Ij slightly below its 
critical current Iq , the rate of switching is therefore very 
sensitive to noise in the current. The first detection of 
asymmetric noise with a J J was reported in Refs. 0,3- 
However, the detector JJ, which was placed in an induc- 
tive environment, had a very large plasma frequency, and 
the dynamics of the junction changed regime as the noise 
intensity increased, from macroscopic quantum tunneling 
(MQT) to retrapping fld\ through thermal activation. 
The measured asymmetry in the escape rates could only 
be compared to an adiabatic model (llj . using empiri- 
cal parameters. For a detector to be of practical use it 
must have a well characterized and a simple enough re- 
sponse, so that information on the noise can be reliably 
extracted. As quantitative theories have been developed 
for a JJ injhe_thcrmal regime placed in a resistive envi- 
ronment jl2l . [isl ll^ , we designed an experiment in this 



framework, allowing for a detailed, quantitative compar- 
ison with theory. 

The principle of our experiment is to add the cur- 
rent noise from a noise source to the DC bias current 
of a J J (see Fig. 1). The dynamics of a J J placed in 
a resistive environment are described by the celebrated 
RCS J model , with the voltage related to the average 
velocity of a fictitious particle placed in a tilted wash- 
board potential. The tilt of the potential is determined 
by the reduced parameter s = Ij/Iq. At s < 1 the po- 
tential presents local minima where the particle can be 
trapped. The voltage is then zero: this corresponds to 
the supercurrent branch. The frequency of small oscilla- 
tions is called the plasma frequency ujp. Johnson-Nyquist 
current noise related to the finite temperature T of the 
environment of the junction is modeled as a fluctuat- 
ing force on the particle, which triggers escape from the 
local minimum ("switching"). When ksT > hiOp /2tt^ 
the switching rate F is given by Kramer's formula [T^ 
r = Aexp ~B2(T) with A ~ ujp/2'K for moderate quality 

factors Q, ujp = ujpQ (l — s^)^^* the plasma frequency in 
the tilted potential, ujpo = y^Io/tpQCj the bare plasma 
frequency determined by the capacitance Cj and critical 
current Iq of the junction, with ipQ = h/2e, and 



B2{T) - (4V2Ioipf,/3kBT) (1 - |s|) 



^3/2 



(1) 



Recently, this result was extended to the situation in 
which an additional delta-correlated noise SlN{t), char- 
acterized by a finite third cumulant ^3 defined by 
{SlN{t)SlN{t')SlN{t")) = S36 {f - t) S {t" - t') and a 
second cumulant (SlN{t)6I]y(t')) = S2S {t' — t) adds to 
current through the J J [ij, [ij, [3 • The effect of higher 
order cumulants is assumed to be weak. The corre- 
sponding fluctuating force leads to a modification of the 
rate: F = Aexp — (i32(Teff) -I- B3) . The second cumulant 
yields an increased effective temperature Teff given by 



2fcBreff/i?|| = 2kBT/R + S2 



(2) 
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Here, R is the parallel combination of all the resistances 
which produce Johnson-Nyquist noise across the junc- 
tion. The resistance R» characterizes the friction acting 
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Vn (mV) Vj (^V) 

FIG. 1: (color online) Top: Experimental setup. Noise from 
a tunnel junction (green double box) biased at DC with 
a current In couples through capacitors Ci — 230 pF and 
C2 = 345 pF to a Josephson junction detector (J J, in orange), 
which is current-biased on its supercurrent branch. The volt- 
age Vj across the junction monitors the switching to the dis- 
sipative state. Capacitor Cj = 12.5 pF lowers the JJ plasma 
frequency to LJpo/2'K ~ 1.5 GHz. Capacitors Cjv = 190 pF 
and Cg = 140 pF shunt the external impedance at ujpo, so 
that the impedance across the JJ is determined only by on- 
chip elements. Resistors Ri = R3 — 215 Q and R2 = 515 
were fabricated with thin Cr films. Bottom left: IV charac- 
teristic of the tunnel junction, linear at this scale, with inverse 
slope Rn = 22.9 k57. Bottom right: IV characteristic of the 
detector JJ, with critical current Jo = 437 nA. We attribute a 
resonance near Vj ~ 120 /iV to a mode of the electromagnetic 
environment of the junction. 

on the fictitious particle and is, in a simple model, given 
by the total resistance across the junction, including both 
R and the resistance Rn of the noise source. This expres- 
sion indicates that the second cumulant of noise from the 
noise source S2 simply adds to the Johnson- Nyquist noise 
of the rest of the circuit. The third cumulant gives rise 
to the additional term 

B3 = -S3 {(po/kBTefff ujIq j{s) (3) 

with j{s) a function of the tilt that depends on the quality 
factor [U] . When reversing the sign of the average cur- 
rent In through the noise source, S2 remains unchanged 
whereas S3 changes sign. Therefore, the departure from 
1 of the rate ratio 

i?r = r( + /Ar)/r(-/Ar) = CXp(2|B3|) (4) 

is a measure of non-symmetric noise (5*3 ^ 0) . 

The experimental setup is shown schematically in 
Fig. 1. As it is well established that current noise through 
a tunnel junction is Poissonian (5*2 — e\lN\ and S3 — 
c^In, with e the electron charge), we use such a device 



(green double box) as a benchmark noise source. The de- 
tector J J (orange crossed box) is coupled to it through ca- 
pacitors Ci and C2 . The finite frequency part 51 n of the 
current through the tunnel junction lN(t) flows through 
the detector J J, owing to the high-pass filter formed by 
i?3, Ci and C2 (3dB point at 5 MHz). The switching 
of the J J current-biased at Ij is signaled by the appear- 
ance of a voltage Vj across it. The low plasma frequency 
of 1.5 GHz guarantees ksT > hujp/2n even at the low- 
est temperature of our experiment (20 mK) In the 
relevant range of frequencies slightly below u!po/2Tr, nu- 
merical simulations of the actual circuit indicate that the 
quality factor of the Josephson oscillations Q is close to 
5, insuring an underdamped dynamics, and no effect of 
retrapping 0, as long as s 4/7rQ ~ 0.25. 

The sample was fabricated on a thermally oxidized 
high resistivity (lO'^ to lO'^fJcm) Si wafer. All on-chip 
resistors are lOnm-thick Cr layers, with 215^l/D sheet 
resistance at 4K, placed between mm-size pads. Capac- 
itors were obtained from parallel aluminum films sepa- 
rated by 29 nm-thick sputtered silicon nitride as an insu- 
lator [ISi] . The tunnel junction and the detector J J were 
fabricated at the same time by shadow evaporation of 
20 nm and 80 nm-thick aluminum films. Their current- 
voltage characteristics are shown in Fig. 1. The tunnel 
junction has an area of 0.09 jiui^ and a tunnel resistance 
Rn = 22.9 kJl. It was biased at voltages larger than twice 
the superconducting gap 2A/e = 0.4 mV (which corre- 
sponds to In = 0.02 fiA), so that it behaves as a normal 
metal junction, with Poissonian noise. The detector JJ, 
with area 1 /^m^, has a supercurrent Jq = 0.437 ^A. It was 
biased in series with a resistor Ri — 215 through a 50 H, 
coaxial line equipped with attenuators. When switch- 
ing occurs at a supercurrent Isw, the voltage across the 
junction jumps to (i?i -|-50rJ)/sw < 2A/e, so that the 
current through it drops to zero and no quasiparticles 
are generated. Moreover, gold electrodes in good contact 
with the Al films were fabricated a few /im away from 
the junctions in order to act as traps for spurious quasi- 
particles [l^ that could be excited by the high frequency 
noise. Apart from the Cr resistors and the Au traps, all 
conductors on the chip are superconducting aluminum 
films. 

The sample was thermally anchored to the mixing 
chamber of a dilution refrigerator. The tunnel junc- 
tion was biased by a fioating voltage supply through two 
1.5 Mfi resistors. The on-chip capacitance Cn = 190 pF 
on the bias line is large enough to maintain the voltage 
across the tunnel junction at V^v = RnIn for all rele- 
vant frequencies. Escape rates of the JJ were measured 
using 2 X 10^ current pulses of duration t = 0.53 fis with 
alternatively positive {+Ij) and negative {—Ij) ampli- 
tude, separated by 9/iS. They were fed through a non- 
polarized capacitor Cin — 200 fiF placed at room temper- 
ature, which prevents DC thermoelectric currents from 
unbalancing the pulses. The switching rates r_|_ and r_ 
for the two signs of /,/ were deduced from the switching 
probability P — 1 — e"'"'^ measured as the fraction of the 
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FIG. 2: (color online) (a) Dependence of B'^^^ = 
(— log(r/A))'^'''' on the current /,/ through the detector JJ, 
for data taken at T = 20 mK and currents through the tunnel 
junction 7jv = 0.02 to 2.11 /jA, by steps of 0.42 ^lA. Linear 
dependence is a signature of the thermal activation regime, 
(b) Effective temperature extracted from the slope of datasets 
as in (a), as a function 7jv, for various temperatures T. Ar- 
rows indicate the datapoints corresponding to the plots in (a). 
Solid lines are linear interpolations. Dashed line are the pre- 
dictions from full theory, taking into account the frequency 
dependence of the admittance Y{ijj) across the JJ and of the 
transfer function aiua) from noise source to JJ detector. 



current pulses which led to a voltage pulse. 

We first demonstrate that the switching of the detector 
junction is well described by the model of thermal activa- 
tion whatever the current in the noise source. Figure 2 (a) 
shows, for various currents In > 0, the s— dependence of 
^2/3 = [-ln(r/A)]^/^ Data fall on straight lines that 
extrapolate to for Ij = Jq, as expected from Eq. ([T]). 
This allows us to extract an effective temperature Teff, 
whose dependence on In is shown in Fig. 2(b), with data 
taken at four different base temperatures T. We do find a 
linear dependance with correct extrapolations a,t In — 
(values slightly above T are attributed to imperfect filter- 
ing), as expected from Eq. ^ with 52 = el. Understand- 
ing the slope quantitatively requires an accurate model 
of the actual circuit at microwave frequencies: the RCSJ 
model assumes that the JJ is simply connected to a ca- 
pacitor, a resistance R\\ and a current source, in parallel. 
In the limit Q ^ R\\, which describes friction, has to be 
replaced with _R||(wp) = l/Re(F(wp)), with Y (cj) the 
total admittance of the circuit across the J J. Microwave 
simulations indicate that R\\ (wp) varies almost linearly 
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FIG. 3: (color online) (a) Rate asymmetry 7?r — 1 as a func- 
tion of current in noise source In for various temperatures T. 
Filled (open) symbols correspond to measurements with pos- 
itive (negative) current through the detector JJ. Solid lines 
are theoretical predictions scaled by a factor 1.5. (b) Depen- 
dence of Rr — 1 on s = Ij/Io, at In = O.SSAfiA, obtained 
with measurement pulses of various lengths. Solid line is the- 
oretical prediction scaled by 1.5. 



from 63 at 1 GHz to 36 at 1.5 GHz, and that a cur- 
rent /Ar(ci^) through the tunnel junction leads to a current 
a(a;)/Ar(ci;) through the detector J J, with a transfer func- 
tion a{uj) varying from 1.1 at 1 GHz to 1.27 at 1.5 GHz. 
Since escape is determined essentially by the noise at Wp, 
we replace 5*2 by a'^{u!p)elN ■ Altogether, the prediction 
Teff ~ r -I- a'^{ujp)R^\{ujp)elN /'^kB is in agreement with 
the data (see dashed lines in Fig. 2; to fit the 20mK-data, 
we used T ~ 72 mK), apart from the slight change in 
slope when varying T which could be attributed to vari- 
ations in the kinetic inductance of the superconducting 
electrodes. 

We now discuss the effect of noise asymmetry. The 
^2/3 piQ^g fQj- opposite signs of the current through 
the noise source are undistinguishable within the symbol 
size, demonstrating that the effect of the second cumu- 
lant S2 is dominant. In the limit cVn ^ ksT, theory 
predicts that the effect of S3 is to shift the curves by 
AIj ~ 0.6Biujpo/Q){S3/S2) = 0.6Be/RnCj ~ 0.2% x/o, 
which is difficult to measure reliably [20]. In our experi- 
ment, we measured directly the asymmetry ratio Rr de- 
fined by Eq. (H]), which varies by several % (see Fig. 3). 
We first set the amplitude Ij of the current pulses at a 
value corresponding to a switching probability P ~ 0.6, 
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for which the statistical precision on the rates is good 
PH . We then measured 100 times and r_, with al- 
ternatively +In and —In through the noise source. This 
allows for two independent measurements of i?r : — 

TU-In)/T+{+In) and i?f = {+In) /T_ {-Im) 
[22| . In Fig. 3(a), we show with full and open symbols 
the corresponding measurements. The rate ratio i?r dif- 
fers from 1, a signature of asymmetric noise, as soon as 
In ^ 0- The statistical uncertainty on i?r is smaller than 
the symbols. Small differences between i?p and i?p, in 
particular around In = 2 /uA, are not understood. As 
for the comparison with theory, a difficulty arises be- 
cause of the frequency dependence of the transfer func- 
tion a(a;), which results in a colored third cumulant at 
the detector S'3(a;i,a;2) = a{uji)a{uj2 — LL!i)a{—LU2)e.'^lN- 
In the following, and in the absence of indication as to 
which frequencies are important, we compare however 
with the only existing theory, which assumes white noise 
{S3 = c^In)- The corresponding predictions, Eqs. (|3|4p 
with i{s) ~ 0.81 (1 — s)^'^^[23j shown as solid lines 
in Fig. 3(a), scaled by an arbitrary factor 1.5. Rr exhibits 
a maximum as a function of In due to the opposite vari- 
ations of 5*3 and T^s with In- For Teff, we used interpola- 
tions between the measured values shown in Fig. 2. When 
scaled up by 1.5, which might be due to frequency depen- 
dent transmission [a{uj)), theory accounts well for the 
experimental data. Feedback corrections due to the de- 
tector, described in [l^, are neglected since R\\/Rn <C 1 
[T^ . Note that there is no feedback associated to the 
series resistance i?3 like in Ref. Q because the current 
noise associated to does not flow through the noise 
source, but through the detector J J. In Fig. 3(b), we also 



compare with theory the s— dependence of i?r- In or- 
der to perform this measurement, we used pulses of vari- 
ous durations (0.53 to 21 /xs), which allows to obtain the 
switching rates at different values of s. For the longest 
pulses, the rate asymmetry is as large as 16%. Here also, 
theory scaled by 1.5 accounts precisely for the data. 

Qualitative agreement between experiment and theory 
gives confidence for the use of the J J as a measuring de- 
vice for 5'3, even if the application to a wider range of 
systems requires some theory for colored noise. A limi- 
tation concerns situations with strong non-linearities in 
the voltage dependence of the cumulants, where feedback 
effects could become sizeable [l^, HI]. For quantitative 
measurements of 53 on other systems, it is not only im- 
portant to tune the plasma frequency of the junction in 
the GHz range as done in this work, but also to improve 
the microwave design, in particular with more compact 
electrodes, so as to avoid frequency dependent factors in 
the analysis. Proposals to access the full counting statis- 
tics with a JJ embedded in more complex circuits Q 
remain to be investigated. 
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